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Abstract: Surface properties of neutron-neutron (T=1) pairing in semi-infinite nuclear 
matter in a hard wall potential are investigated in BCS approximation using the 
Gogny force. Surface enhancement of the gap function, pairing tensor and 
correlation energy density is put into evidence. 

The possibility of producing exotic nuclei has spurred an intense revival in 
nuclear structure research. Among many interesting properties neutron-neutron (T=1) 
pairing is certainly one of the main issues in nuclei with a large neutron skin on the 
one hand or neutron-proton pairing (T=0) in heavier N=Z nuclei on the other hand. In 
this letter we want to concentrate on n-n pairing postponing the n-p pairing to future 
work. 

In finite nuclei pairing properties are strongly fluctuating due to shell effects 
and the amplitudes of the fluctuations of the gap are of the same order as its mean 
value [1]. Even Strutinsky averaging does not help very much to obtain for example a 
non fluctuating pair density as a function of the nuclear radius. Another way, which is 
often used to obtain smooth properties free of shell effects, is to consider semi-infinite 
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nuclear matter. This is for instance tine case for tine study of surface properties and 
that is also here the way we will follow to investigate superfluid properties in the 
nuclear surface. Since fully self-consistent HFB calculation for semi-infinite nuclear 
matter are extremely difficult compared to the non superfluid case [2], we will use in 
this exploratory work a rather simplified model which, however, should allow to reveal 
the essential features. To this purpose we embed nuclear matter in a semi-infinite 
hard wall potential and solve the pairing problem in BCS approximation with the 
Gogny force. The finite range Gogny force (D1S) [3] gives realistic pairing properties 
in finite nuclei [3]. In the S=0 T=1 channel it is density independent and yields in 
homogeneous matter practically identical gap values as those calculated with the 
Paris force [1,4]. Therefore, the Gogny force acts for pairing properties like a free 
force, a conclusion that also has been found in ref [5]. In the past it often has been 
invoked on indirect grounds that nuclear pairing acts predominantly in the nuclear 
surface. However, never any extended investigation of the r-dependence of pairing 
quantities such as the gap, the anomalous or pair density, or the correlation energy 
density has been undertaken using a realistic pairing force for the semi-infinite 
nuclear matter problem. Only in the context of inhomogeneous nuclear (neutron) 
matter in the inner crust of neutron stars some limited studies on the r-dependence of 
pairing properties have recently been undertaken [6]. The conditions are, however, 
quite different from the ones of real nuclei. 

Let us therefore define the non-local gap in r-space (we suppress the spin 
indices), 

A(r,r') = 2;^,(r)^_,(r')A(k) (1) 

k 
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where k is the momentum and A(k) is the gap which in BCS approximation is 
diagonal and ^^(r) are the wave functions corresponding to the one sided infinite 
wall potential (z is the coordinate normal to the wall), 



(p^{r) = ©{z)sm{k^z)e-'^'^ (2) 
With (2) we obtain for(1), 

A{ry) = -^&{z)@{z{\dpSoos{p^{z - z))-cos{p^{z + z))} 

(3) 



where Jg 's the Bessel function of zero order. A useful quantity is the Wigner 
transform of the gap, 

A(R,k) = jjsA(R,s>-*^ (4) 

where, R = r + r'/2and s = r-r' are center of mass and relative coordinates. From (4) 
we directly obtain. 



^ - ^^^-'^^ (5) 

It is interesting to recover from (5) the value of the gap in the bulk in considering R^- 
values far from the wall, 



lim A(i?2,k)= \ dp^S{p^ -kjA{p^,kJ-2S{kJ lim \dp^cos{2p^R^)\{p^,kJ 

„ (6) 
= A(k^,k^)-ld{k^) lim \dp^cos(lp^R^)A[p^,k^) 

Due to the strong oscillations of lim cosilp^R^) the second term of (6) will vanish so 
that we find the important simplifying result, 



lim A(i?^,k) = A^(k) = A(k) (7) 



where the index B stands for "bulk". Therefore we can replace A(k) by A.g{k) in all 

preceding formulas. This is very helpful, since we have shown previously [7] that the 
following analytical approximation to the solution of the gap equation in the bulk is 
quite accurate for the Gogny force over the whole range of k-values,* 

A.(^) = ^^A,(^,) (8) 
y(kp,kj,) 

where v(p,k) are the matrix elements (p,-p|v|k,-k) averaged over the angle 
p.k/pk which is the s-wave part of the force in momentum space and reads. 



1 1 -M^^-^ .J ,pk^ 



v(p,fc) = Fz, — — e ' " sinh 



(9) 



*) In [7] only t^^k,) had been tested (see eq. (10) below). We, however, checked that (8) is well 



verified also for k i^kp. 
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where z, =n'^idl{yv^ -B^-H^ +Mjwith the values of the parameters given in [3]. 
The bulk gap at the Fermi level is given by [7], 



(10) 



with A^(0) = ^ 



m 



being the level density at ^^(no spin, no isospin factors), 



m* =m*{kp) the effective mass corresponding to the Gogny force and, 
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(11) 



where, 



g{x) = -N{0)x 



/(^xkp , k p)w{k p , xk p ) 
\{kp,kp) 



(12) 



Another interesting quantity is the pairing tensor. In momentum space it is given by, 
A(^) 



4k) = 



IE,. 



(13) 



where. 



^ e ki 



2m 2m 

y 

is the usual quasiparticle energy 



(14) 
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With Ag{k) from (8) it is tlien easy to also calculate the Wigner transform /f(/?^,k) of 
the pairing tensor. Integrating A(^^,k)and }c{R^,k) over k yields the corresponding 
quantities in r-space for r - r' = 0, that is the local part of the gap and pairing tensor. 
One obtains, 



f \ 
2 
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sin,... 



(15) 



2 ^ll] 1 



\kdk- 



1^ r 
e ^ sinh 



2 kk p 



2E,. 



sin(2iti?J 
2kR, 



(16) 



with, 



D = 



( 2 ^ A,(^,) 1 



(17) 



In Fig. 1 we show these two quantities together with the density. 



sin(2 pR_) 



i{k) 



(18) 
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where, 




(19) 



are the BCS occupation numbers. We see that A(i?J has a strong peak very close to 
the wall but is otherwise rather structureless. On the other hand in J rather 
pronounced Friedel oscillations are apparent but no dramatic surface enhancement 
occurs. One notices however that with respect to p{R,) the anomal density >c{rJ is 
pushed closer to the wall. However, A(i?J and }c{rJ contain only limited information. 

The quantity, which really measures the importance of pairing correlations, is the 
correlation energy density. 



In Fig. 2 we show £^„„{rJ as a function of R^. We see that ^^^^^ follows a similar 
pattern as k{r^): the last bump is relatively more important than in p{Rj and ^^<„,is 
closer to the wall than p{Rj. This is probably a combined effect of Friedel oscillations 
and the fact that Ag{kp) increases when kp, i.e. the density decreases. Another 
quantity, which is sensible to be studied in the present context, is the gap function. 



This stems from the fact that if(^,,p) is very much peaked for \p\ = kp{Rj where 



approaching the surface it remains true approximately. This can be seen from the 




d 



(20) 



A(R^,\p\ = kp{Rj)^Ap{Rj 



(21) 



kp{R Joe p^^^{rJ is the local Fermi momentum. This is true in the bulk. When 
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product A{R^,\p\ = kp{R^)).K:{R^) which approximates s^^rri^) rather well (dotted curve 
in Fig. 2). The quantity A^(i? J is shown in Fig. 3. One again notices a quite strong 
surface enhancement. 

All in all one finds that the surface enhancement of pairing with respect to its bulk 
value is relatively moderate, not exceeding 20% to 30%. These values may increase 
somewhat for a more realistic Saxon-Wood potential. For an extremely smooth 
potential the local density approximation (LDA) should be valid. The latter yields 
considerably more surface enhancement [1]. Realistic potentials are smoother than 
the hard wall potential used here. However, they are by far not smooth enough to 
make LDA an exact theory. All that one can hope is that average quantities of finite 
nuclei are well reproduced. This is what happens, since e.g. correlation energy and 
level density enhancement compare well with quantal calculations [1 ,8]. On the other 
hand one should realize that the gap values found in the present case of half in finite 
matter are smaller by an important factor compared to those of finite nuclei. Whether 
this enhancement, which comes essentially from the strong size dependence of the 
matrix element of the pairing force on the Fermi surface, is a volume or a surface 
effect remains to be seen by future studies. 

The situation may also change for exotic nuclei close to the drip line with a neutron 
skin. First of all one knows [9] that in this case HFB rather than BCS must be used. 
Second our potential certainly can not describe the situation of drip line nuclei. 
Nonetheless we do not think that even in this case the situation changes dramatically 
with respect to our present study. It seems very unlikely that for example the 
correlation energy in the surface can attain values 2 or 3 times bigger than in the 
bulk. 



A further interesting fact, wliicli can be revealed from our study, is tliat tine 
momentum dependence of pairing quantities becomes distorted close to the surface. 
Such features have already been revealed for the normal density [10] but here it 
repeats itself for the anomal density. Indeed the quadrupole moment 
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+A',{k) 



(22) 



7m [k) 2m* {kp) 

shows quite strong (positive) deviations from zero close to the wall (Fig. 4). It means 
that the anomal momentum distribution has a prolate deformation with the long axis 
perpendicular to the wall. The S'{r^) term is an artifact of the hard wall potential. In 
any smooth potential the J' would be smeared out filling in the negative dip in Qjjust 
before the location of the wall. The prolate momentum distortion is physically quite 
understandable, since the wave function of the Cooper pair has more freedom 
(space) parallel to the surface than perpendicular to it (Heisenberg's uncertainty 
principle). The anisotropic momentum distribution of the Cooper pair can influence 
pair transfer reactions and, in condensed matter physics, tunneling rates in 
Josephson junctions. 

Having the surface behavior of nuclear pairing at hand one can ask the 
question how pairing influences the surface tension cr. Since the pairing is relatively 
more important in the surface, one expects that this extra binding leads to a reduction 
of the surface energy in spite of the fact that the slightly enlarged surface profile due 
to A leads to an increase of a in the normal part of the energy. This effect is however 
extremely small and one indeed finds that the surface tension is lowered due to 



pairing. However for n-n pairing tliis is only in tine 1% range. Since n-p pairing is 
supposed to be considerably stronger [11] one can assume that the nuclear surface 
tension is lowered by maximally 10% for heavier N ~ Z nuclei. 

In conclusion we have shown that surface enhancement of nuclear pairing is 
important but not dominant. Our studies are based on the Gogny D1S force, which 
gives in nuclear matter quite similar gap values as the Paris force. Our simplified 
model of a half-infinite hard wall potential should not invalidate our general 
conclusions. The local density or Thomas Fermi approximation seems to indicate that 
the surface behavior is considerably more important than the one found here. One 
should, however, remember that these semi-classical solutions should be treated as 
distributions and that expectation values like e.g. the correlation energies or level 
densities, in comparison with their quantal counterparts, are surprisingly well 
reproduced [1 ,8]. 

NOTE: While this work was being completed we became aware of a preprint by M. 
Baldo et al. (nucl-th 9812018) containing qualitatively similar results on the surface 
behavior of nuclear pairing. 

Acknowledgments: we thank M. Baldo, E. Saperstein, X. Vihas for discussions. 
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Figure captions 



Figure 1 : the local part of the gap A(i?J, the pairing tensor }c{rJ and the 

density p{Rj, normalized to their bulk (infinite nuclear matter) values, 
respectively -.84 MeV, .012 MeV and .1632 fm'^as functions oiR^. 



Figure 2 : the correlation energy density e^^^^{Rj , dotted curve, and the product 

A(R^,\p\ = kp{R^)).K{R^) normalized to the bulk values (^^^^^(oo) = -008 
MeV) as functions oiR^ . 

Figure 3 : the gap function a{r^ , |p| = hp {r^ )) = Ap{Rj and the 

density p{rJ, normalized to their bulk values (Ap(+oo) = o.60Mey ). 



Figure 4: the quadrupole moment Qjl^Jof the Wigner transform if(^,,k) of 
the pairing tensor in semi-infinite nuclear matter 
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